hep-ph/9610370| SLAC-PUB-7316 

Saclay/SPhT-T96/lll 
UCLA/96/TEP/33 
October, 1996 



One-Loop Amplitudes for e+e — > qqQQ 



Zvi Bern" 

Department of Physics and Astronomy 
University of California, Los Angeles 
Los Angeles, CA 90095 
bern@physics.ucla.edu 

Lance Dixon* 
Stanford Linear Accelerator Center 
Stanford University 
Stanford, CA 94309 
lance@slac . Stanford. edu 

David A. Kosower and Stefan Weinzierl 
Service de Physique Theorique* 
Centre d'Etudes de Saclay 
F-91191 Gif-sur-Yvette cedex, France 
kosower@spht . saclay . cea.fr 
stef anwOspht . saclay . cea.fr 



Abstract 

We present the one-loop helicity amplitudes for processes involving a vector boson V (V = 
W,Z, or 7*) and four massless quarks, — ► VqqQQ, where V couples to a massless lepton pair. 
These amplitudes are required for next-to- leading order 0{oP s ) numerical programs for four-jet pro- 
duction at e + e~ colliders, for W, Z or Drell-Yan production in association with two jets at hadron 
colliders, and for three-jet production in deeply inelastic scattering experiments. We obtained the 
amplitudes presented here by using their analytic properties to constrain their form. 
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1. Introduction 



Electron-positron annihilation provides a clean experimental laboratory for studying jet prop- 
erties. Leading-order predictions for the production of up to five jets have been available for quite 
some time [1,2,3,4], but the reduction of theoretical uncertainties requires next-to-leading-order 
(NLO) QCD corrections. The NLO matrix elements for three-jet production and other 0(a s ) ob- 
servables are also known [3] , and numerical programs implementing these corrections [5] have been 
widely used to extract a precise value of a s from hadronic event shapes at the Z pole [6]. 

Next-to-leading order corrections for more complicated processes are required, however, if 
we wish to use QCD in probing for new physics in other standard model processes. In e + e~ 
annihilation, for example, four-jet production is the lowest-order process in which the quark and 
gluon color charges can be measured independently. Four-jet production is thus sensitive to the 
presence of light colored fermions such as gluinos [7]. At LEP 2 the process e + e~ — > (7*, Z) — ► 4 jets 
is a background to threshold production of W pairs, when both Ws decay hadronically. The one- 
loop matrix elements required for an NLO study of four-jet production are also needed for a 
next-to-next-to-leading (NNLO) study of three-jet production at the Z pole. Such a study (which 
awaits the computation of certain two-loop matrix elements as well) would be desirable in order to 
reduce the theoretical uncertainties in determining a s via this process. 

In this paper we present analytic formulas for the one-loop helicity amplitudes for electron- 
positron annihilation into four-quarks, e + e~ — ► (7*, Z) — > qqQQ. Together with the leading-in-color 
matrix elements for production of two quarks and two gluons, e + e~ — > qqgg [8,9], the leading-color 
parts of these amplitudes have already been incorporated into an NLO program for e + e~ — ► 4 jets 
[10]. The same amplitudes presented here may also be used in computations of W or Z + 2 jet 
production at hadron colliders and three-jet production in deeply inelastic scattering. Glover and 
Miller [11] have also recently reported on a calculation of the squared one-loop matrix elements for 
e + e~ — > 7* — > qqQQ, summed over helicities and expressed in terms of Lorentz scalar products of 
the quark four-momenta, rather than the spinor products that we employ. A comparison of their 
results with those presented in this paper will be useful. 

Recent years have seen a number of technical advances in the computation of one-loop am- 
plitudes, which the authors have surveyed in a recent review article [12]. These advances have 
made possible the calculation of all one-loop five-parton processes [13,14,15], as well as of a number 
of infinite sequences of one-loop amplitudes [16,17,18,19]. The general strategy employed in this 
paper (and in a subsequent paper on e + e~ — > qqgg [8,9]) is to obtain amplitudes from their analytic 
structure. In particular, we use the constraints of unitarity [20,18,19,21] and factorization [16,22], 
as summarized in ref. [12]. This approach leads to relatively compact expressions, as compared 
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with those obtained from a traditional, diagrammatic computation. In this approach, amplitudes 
obtained previously are recycled to obtain new amplitudes; manifest gauge invariance is therefore 
maintained. In a Feynman diagram approach each diagram alone is not gauge invariant, and are 
often individually much more complicated than a final sum over diagrams. 

We use helicity methods [23,24] since they lead to relatively compact expressions for the am- 
plitudes, and retain all spin information. We also make use of color decompositions [25] to help 
simplify the analytic structures that must be computed. As a check, we have verified numerically 
that the amplitudes presented in this paper agree with a direct Feynman diagram calculation we 
performed. 

The paper is organized as follows. In section 2, we briefly describe helicity methods and color 
decompositions. We give the amplitudes, together with a brief description of the calculational 
methods, in section 3; we describe the contribution proportional to the axial vector coupling of the 
Z to the t, b quark isodoublet in subsection 3.4. A summary is included in section 4. We collect 
descriptions of various integral functions appearing in the amplitudes in an appendix. 



2. Basic Tools 

In this section, we briefly review two of the basic tools useful for expressing amplitudes in 
a compact form: the spinor helicity method and color decompositions. The reader is referred to 
review articles [26] and references therein for further details. 

2.1 Spinor Helicity 

In explicit calculations it is usually convenient to use a helicity basis, where all quantities 
are rewritten in terms of Weyl spinors |/c ± ). Although there are no external gluons in the final 
e + e~ — ► qqQQ helicity amplitudes, they appear as intermediate states in various unitarity cuts 
and factorization limits that are used to construct the amplitudes. We made use of the gluon 
polarization vectors of Xu, Zhang and Chang [23,24], 

mV ' V2{qk) ^ ' V2[kq] K ' 

where k is the gluon momentum and q is an arbitrary null 'reference momentum' which drops out of 
final gauge-invariant amplitudes. The plus and minus labels on the polarization vectors refer to the 
gluon helicities. Our (crossing-symmetric) convention takes all particles to be outgoing, and labels 
the helicity and particle vs. antiparticle assignment accordingly. (That is, we write the amplitudes 
for the process — > VqqQQ.) For incoming (negative energy) momenta the helicity and particle 



vs. antiparticle assignment are reversed. It is convenient to define 

(ij) = (kr\k+), [ij] = (kt\kr), 

(i\l\j) = (kr\^\kj), (i\(i + m )\j) = ( k -\^ l + it m )\kr), (2.2) 

+ m)(n + r) • • • = (k~ | (ft + f m ){f n + %) ■ ■ ■ | fc+) , 

for the spinor products, which is the notation we shall use to quote the results. Here all the 
momenta ki are massless. The spinor inner products (ij), [ij] are antisymmetric and satisfy 
(ij) [ji] = • kj = Sij. 

To maximize the benefit obtained from the spinor helicity formalism for loop amplitudes we 
must choose a compatible regularization scheme. In conventional dimensional regularization [27], 
the polarization vectors are (4 — 2e)-dimensional, which is incompatible with the spinor helicity 
method's use of four-dimensional polarizations. To avoid this problem, we modify the regularization 
scheme so all helicity states are four-dimensional and only the loop momentum is continued to 
(4 — 2e) dimensions. This is the four-dimensional- helicity (FDH) scheme [28], which has been shown 
to be equivalent at one-loop [29] to an appropriate helicity formulation of Siegel's dimensional- 
reduction scheme (DR) [30]. The conversion between schemes has been given in ref. [29]; there is 
no loss of generality in choosing the FDH/DR scheme. 

2.2 Color Decomposition 

In this section we describe a color decomposition [25] of the one- loop amplitude for e + e~ — ► 
qqQQ, in terms of group-theoretic factors (color structures) multiplied by kinematic functions called 
partial amplitudes. Because of the crossing symmetry of the spinor products [24] and of the integral 
functions (discussed in the appendix), these partial amplitudes may also be used to obtain the one- 
loop contributions to Z or W + 2 jet production at hadron colliders or three jet production in 
deeply inelastic scattering. (For the case of the W, only the coupling constants need be changed 
in the formulae given below.) 

The partial amplitudes are defined to be the coefficients of the various color structures. Con- 
sider the amplitude Ae(l q , 2q, 3q, 4^; 5e, 6 e ). At tree-level its decomposition is 



X ree (l g ,2 ,3 Q ,4,) = 2eV 



(-Q 9 + ^U^( 5 56))< ee (l„2Q,3 Q ,4 g -) 
+ (~Q Q + vl R vl R V z (s 56 ))Al^(3 Q ,4 <1 ,l q ,2 Q ) 



(2.3) 



where we have suppressed the 5, 6 labels of the electron pair, e is the QED coupling, g the QCD 
coupling, Q q (Q®) is the charge of quark q (quark Q) in units of e, and the left- and right-handed 



couplings are 



,,e 



-l + 2sin 2 6V „ 2sm 2 6 w 



Vl sm.20 w ' Vr s\n26 w ' ^ ^ 

q ±1 — 2Q q sin 2 0vk .. 2Q«sin 2 w 

Vl ~ sin 26 w ' Vr ~ sin 29 w ' 

where 6w is the Weinberg angle. Eq. (2.3) contains the ratio of the Z and photon propagators, 

= m| + (2 ' 5) 

where and are the mass and width of the Z. We have given the decomposition for a general 
SU(N C ) gauge group; N c = 3 for QCD. The two signs in v q L correspond to up (+) and down (— ) 
type quarks. 

The subscripts L and R refer to whether the particle to which the Z couples is left- or right- 
handed. That is, v q R is to be used for the configuration where the quark (leg 1) has plus helicity and 
the anti-quark (leg 4) has minus helicity, which we denote by the shorthand (1^~,4^~). Similarly, v q L 
corresponds to the configuration (l~,4j~). Because the electron and positron are incoming in e + e~ 
annihilation, our outgoing-momenta notation reverses their helicities and particle vs. anti-particle 
assignment. Thus, v e R corresponds to the helicity configuration (5g , 6+) whereas v e L corresponds to 
the configuration (5g~ , 6~ ) . 

We have defined the tree and one-loop partial amplitudes ylg ree and A 6;i to include a photon 
propagator. The ratio Vz{sb&) appearing in eq. (2.3) then replaces the photon propagator with 
a Z propagator. This form of the amplitude is convenient for checking that amplitudes properly 
reduce to lower-point amplitudes when the e + and e~ momenta are taken to be collinear. 

At one loop there are three partial amplitudes, 

"4-6 P (lg> 2q> 3q, Aq) = 



2eV 



{~Q q + v e L , R vl R Pz(s 5e )) [N c ^ 6g A 6;1 (l q , 2 Q ,3 Q ,4 q ) + ^ 5£ A 6;2 (l q , 2 Q ,3 Q ,4 q ) 
+ (~Q Q + vI,rvIr V z (s 56 )) [N c 5^ 6£ A 6;1 (3 , 4„ 1„ 2 Q ) + 8£ 5^ A 6;2 (3 Q , 4„ 1„ 2 Q ) 
+ slfc*<«»> ('* C " WS C)^3(1 9 ,2 ,3 Q ,4 9 -) 



(2.6) 

Note the additional factor of N c in the color-tensor coefficients of A e . 1:2 compared to the correspond- 
ing tensors in the tree-level amplitude. The A e . 3 term arises from a fermion triangle graph [31]. 
It violates the axial symmetry, and is proportional to the axial coupling of the Z to the top and 
bottom quark isodoublet as well as to the top-bottom mass splitting (see section 3.4). The cor- 
responding charges vanish for the photon and the W boson, and this term does not contribute to 
amplitudes for the latter bosons. 



Helicity conservation for massless quarks and leptons ensures that there are only 2 3 = 8 possible 
helicity configurations for e+e~ — > qqQQ, namely two choices for each fermion line. Because the 
electron line couples through the current (5 ± |7 M |6 ± ) = (6+|7 M |5+), it is trivial to reverse its helicity 
simply by exchanging 5 and 6 in the partial amplitudes Ag ree and Ag-i, and exchanging v e R <-> v e L in 
the prefactors in eqs. (2.3) and (2.6) as discussed above. We can use parity to reverse all helicities 
simultaneously in the partial amplitudes A§ ee , A^i and A &;2 , by complex conjugating all spinor 
products ((ij) <-> [j i])- The axial vector contribution Aq-s is also complex conjugated, but acquires 
an additional overall minus sign from the 75 in the loop. Thus we are left with just two independent 
helicity configurations, which we may take to be As(lq , 2q, 3+, 4^~ ; 5g , 6+) (we shall again suppress 
the 5, 6 labels below). 

The partial amplitudes can be further expressed in terms of 'primitive amplitudes' [15]. The 
primitive amplitudes are gauge-invariant classes of color-stripped amplitudes from which we can 
build the partial amplitudes. (In ref. [15] the primitive amplitudes were defined to have a fixed 
ordering of the external legs, but here we extend the notion to mean gauge-invariant color-stripped 
building-blocks for amplitudes.) It turns out that some signs in the reduction of partial amplitudes 
to primitive amplitudes depend on which of the two helicity configurations one is considering, so 
for clarity we shall explicitly list the two cases separately. Although color decompositions do not 
depend on the helicity choices, these sign differences appear because we have used the symmetries 
of the primitive amplitudes to reduce the number of independent ones required. 

The formulae for A & -i{l^ , 2~t, 3~,4^~) in terms of the primitive amplitudes are 

2j,3 Q ,4r) = ^+(1,2,3,4) + ^X^' ++ (l,2,3,4) - ^ ++(l, 2, 3, 4) 
" lk 2 > 3, 4) + A+~ (1, 3, 2, 4)) + -L A f(2, 3, 1, 4) , 

c 

^e i2 (l+ 2j,3-,4r) = ^+-(1,3,2,4) - 2^1^. ++(1,2,3, 4) + ++(1, 2, 3, 4) (2.7) 

+ -L (A+- (1, 3, 2, 4) + A++(l, 2, 3, 4)) - -L Af(2, 3, 1, 4) , 

c c 

^ 6;3 (l+,2j,3-,4 g -) = ^ x (l,4,2,3), 

where n s is the number of scalars* (n s = in QCD) and n/ is the number of Dirac fermions. 
Within the context of supersymmetry decompositions [13,32,12] it is natural to divide the fermion 
contribution into s and / pieces as we have done here. 

* As in refs. [13,15], each scalar here contains four states (to match the four states of Dirac fermions) so that n s 
must be divided by two for comparisons to conventional normalizations of scalars. 
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For the other independent helicity configuration A 6;i (lg , 2~, 3+ 4^ ), we have, 



A 6;1 (l+,2^,3j,4r) = A+-(l,2,3,4) + \f-A% +"(1, 2, 3, 4) - ^ +-(1,2,3, 4) 

- ^(^+"(1,2,3,4) + A++(l,3,2,4)) - -^(3,2,1,4) , 

c c 

A 6;2 (l+2^,3j,4r) = A 6 ++(1,3,2, 4)- ^^.+-(1,2,3,4) + ^4+"(l, 2 ,3,4) (2.8) 

+ ^ K + (1, 3, 2, 4) + A+- (1,2,3, 4)) + ^4(3, 2,1,4), 

A 6;3 (l+2^3+ 4r) = -Af(l,4,3,2). 



Note that and Ag x appear with a different permutation of the arguments, and the opposite sign, 
compared with the previous helicity structure. The signs in the above expressions also depend on 
the relative phase conventions of the two helicity amplitudes; we have chosen a convention where 
Ag ree '+ _ (1, 3, 2, 4) = — Ag ree '++(1, 2, 3, 4). Of course, consistency between tree and loop amplitudes 
must be maintained. 

A representative diagram for the primitive amplitudes ^4g' +± (l, 2, 3, 4) and ^4g' +± (l, 2, 3,4), 
corresponding to those terms in the partial amplitudes proportional to the number of fermions 
rif, is depicted in fig. 1. There are only two such diagrams. The 'parent' pentagon diagram for 
the leading-color helicity amplitude Aq (1,2, 3,4) is depicted in fig. 2. (By a 'parent' diagram 
we mean a diagram from which all other diagrams in the set can be obtained via a continuous 
'pinching' process, in which two lines attached to the loop are brought together to a four-point 
interaction — if such an interaction exists — or further pulled out from the loop, and left as the 
branches of a tree attached to the loop.) The parent box and triangle diagrams for the subleading- 
color primitive amplitude Af(l,2, 3, 4) are depicted in fig. 3. (We combine the triangle-parent 
contributions with the box-parent contributions because the former are so simple, and the two 
have an identical color structure.) The fermion-loop triangle diagram, proportional to the axial 
coupling of the Z to quarks, is shown in fig. 4. We neglect the u, d, s, c quark masses, and so only 
the t, b quark pair survives an isodoublet cancellation in the loop, thanks to its large mass splitting. 
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7,Z 




Figure 3. The diagrams for the subleading-color component A§ can be divided into (a) 
those with a parent box diagram), and (b) those with a parent triangle diagram. Note that the 
labelling of legs is different from that in the leading-color diagrams. 



+ ({1.2}** {3.4}) 



Figure 4. Triangle diagram for the amplitude component ^4g x , proportional to the axial 
coupling of the Z to quarks. 

Formulas (2.3) and (2.6) apply to the case of unequal quark flavors, q ^ Q. The equal flavor 
amplitude may be obtained from the unequal-flavor formula by subtracting the same formula with 
q and Q exchanged (or q and Q exchanged, but not both), and then setting Q = q in all the 
coupling constant prefactors. (Equal-flavor cross sections also require an identical-particle factor 
in the phase-space measure, which is (^) 2 for e + e~ — ► qqqq.) 

The virtual part of the next-to-leading order correction to the parton-level cross-section is 
given by 



It is a straightforward exercise (which we leave to the reader) to evaluate this color sum in terms 
of partial amplitudes. 



In constructing the amplitudes presented here, we use two basic analytic properties: that their 
imaginary (absorptive) parts be determined from the Cutkosky rules [20], and that they factorize 
on particle poles. These analytic properties of amplitudes have, of course, played an important role 
in field theory for many decades; we make use of recent developments which allow us to to obtain 
complete amplitudes with no subtractions. In order to maximize the efficiency of the computation 
it is useful to perform the computation in the manner described in ref. [12]. Due to the complexity 
of the kinematics for the processes presented here, further techniques are required to minimize the 
appearance of undesired spurious poles; these will be discussed in ref. [9]. 




(2.9) 



colors 



3. The Amplitudes 
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The primitive amplitudes A^ 6 and A§ are proportional to tree amplitudes and are given by 



4 +± (l, 2, 3, 4) = cr4 ee (l+, 2± 3j, 4,; 5 e " , 6+) 
4 +± (l, 2, 3, 4) = c r 4 ree (l^ 2|, 3J, 4,; 5 e " , 6+) 



1 

3e 



-S23 



"«23 



+ 2 



(3-1) 



where ^4g ree is given below, 



1 



cr 



r(i + e )r 2 (i-6) 



(3.2) 



(4vr) 2 - e r(l - 2e) 
and D = 4 - 2e. 

The contribution Aq x is finite (see below). It is convenient to decompose the remaining prim- 
itive amplitudes into divergent (V) and finite (F) pieces, 



^6 



cr 



where 

a = {++, +-, si} 

labels the primitive amplitude under discussion. The quantities A 6 
tree partial amplitudes appearing in eq. (2.3), 



tree, + ± 



(3.3) 



(3.4) 



coincide with the true 



< ee (l+, 2±, 3j, 4,; 5 e ^, 6+) = A?*> +± (1, 2, 3, 4) . (3.5) 

(We continue to omit the arguments of the primitive amplitudes corresponding to the lepton pair.) 
We express the amplitudes in terms of the spinor products (2.2), the Lorentz products Sij — 
(ki + kj) 2 and tij m = (ki + kj + k m ) 2 , and the following combinations of kinematic invariants, 



<5l2 — Si2 — «34 — «56 



C>34 — «34 — S56 — S\2 , 



^56 — «56 — «12 — S34 , 



A 3 = s 12 + S34 + s 56 - 2si 2 s 34 - 2s 34 s 56 - 2s 56 Si2 



(3.6) 



The latter quantity is the negative of the Gram determinant associated with the set of massive 
momenta {k\ + k 2 , k% + k±, k§ + ke}. 

The amplitudes we present are bare ones, i.e., no ultra-violet subtraction has been performed. 
To obtain the renormalized amplitudes in an MS-type subtraction scheme, one should subtract the 
quantity 



c r N c g 2 



1/11 2 rif 1 n s 



jtree 
"™6 ' 



(3.7) 



,3 3 N c 3 N c 
from the the amplitude (2.6). 

We quote the results in the FDH scheme [30,28], but these may easily be converted to the 
't Hooft-Veltman scheme; to do so one would add the quantity 



-crN c g 2 [^-j^)A. 
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(tree 



(3.8) 



to the amplitude (2.6) and change the coupling constant from the non-standard to the standard 
a-jgiT. The conversion between the various schemes is discussed in refs. [29]. 

3.1 The Helicity Configuration q + Q + Q~q~ 

We first give the primitive amplitude ^4 ++ (l, 2, 3, 4) which contributes to the leading color 
part of Afy\{lq , 2~t, 3~,4^"; 5jT, 6+). This amplitude is odd under a 'flip symmetry', which is the 
combined operation of a permutation and spinor-product complex conjugation: 

flip: 1 <-> 4 , 2^3, 5^6, (a b) <-> [a b] , (a - I b |c~) <-> (a + 1 6 |c + ) = (c~ | 6 |a~) . 

(3.9) 

The corresponding tree amplitude for this helicity configuration is 



4 ree < ++ (1, 2, 3, 4) = i [12]{5A) ^ 1 + 2 ^ + i < 3 4 > t 6 ^ < 5 ' ( 3 + 4 )' 2 ) 



S23 "856*123 



■S23 s 56*234 



(3.10) 



We have 



F ++ (l,2,3,4) = (- < 3 K 1+2 >I 6 > 2 1 12 1 2 < 45 >' 2 



(23) [56] ti23<l|(2 + 3)|4} [23] (56) i m (4|(2 + 3)|1) 

x 



(3|(1 + 2)|6) 
[56] (1|(2 + 3)|4) 



/ -S12 -S 2 3 \ . ~2mh 

hs_i — , — ; + Ls.! (s 34 , t 123 ; Sl2, S 56 ) 

v_ C123 -1123 7 

(l|(2 + 3)|6)[12]Lo(fS) (3|4|6)Lo^ 



S56 
23 



2(2 3) [5 6] ti23<l|(2 + 3)|4> 



il23 *123 (2 3) ti23 

((3|2|1)(1|(2+3)|6))" L 



,2 Ll V -S23 



Li 



+ (3|4|6) 2 t? 23 ^ 



( -«56 ) 
V-*123 J 



'23 



123 



-flip, 

where 'flip' is to be applied to all the preceding terms in F ++ . 
3.2 The Helicity Configuration q + Q~Q + q- 



(3.12) 



We now give the result for A§~(1, 2, 3, 4), which contributes to the leading color part of the 
partial amplitude A 6; i(l^, 2q, 3+ 4^~; , 6+). This amplitude is odd under the same flip symme- 
try (3.9) as Aq + . The tree amplitude is 



Af^ +-(1, 2, 3,4) = -i I 13 ] (54) (2|(1 + 3)|6) _ . (24) [6 1] (5|(2 + 4)|3) 



S23'S56*123 

Note that 4 ree ' +"(1,2,3,4) = -A^ ++ (1, 3, 2, 4). 
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523^56*234 



(3.13) 



For the one-loop contributions we have 

,2 \ e / ,,2 \ e \ o / ,,2 



Note that V +_ (l,2,3,4) = V ++ (l,2,3,4). The finite part is 



^-(1,2,3,4) = (- rno1 ^ + ^ ^ 2)2(3K lty. . , ) !■> " S23 



[2 3] <56> ti23<4|(2 + 3)|1> (2 3) [5 6] t 123 (1 3) 2 (1|(2 + 3)|4) > V-W -£ 123 

[13] 2 (45) 2 (3|(1 + 2)|6) 2 (2|(1 + 3)|4) 2 \ ~2 mfe 

1 [23](56)t 123 (4|(2 + 3)|l) + (23)[56]t 123 (l|(2 + 3)|4)(3|(l + 2)|4)2; - 1 1 34 ' 123 ' 12 ' 



+ 



l(si2 -a 3 4) (51) <5|(3 + 4)|2)(1|(2 + 4)|3) 1 (s 12 - s 34 ) [2 3] (5 1) (5 2) 



2 (56)(1|(2 + 3)|4)(1|(3 + 4)|2) 2 2 (56) (1|(2 + 3)|4)(1|(3 + 4)|2) 

1 (*i34 -*234)[34](51)(5 4) [12] (51) (5 2) (1|(2 + 4)|3) 

2(56)(1|(2 + 3)|4)(1|(3 + 4)|2) (56) (1|(2 + 3)|4)(1|(3 + 4)|2) 

1 (51)[6 2](1|(2 + 4)|3)(* 23 4 - 2S34) 1 (5|(1 + 2)|6)(1|(2 + 4)|3) 
+ 2 " (1|(2 + 3)|4)(1|(3 + 4)|2) 2 ~ + 2 (1|(2 + 3)|4)(1|(3 + 4)|2) 

(l|(2 + 4)|3)( Sl2 5 12 ((5|l|6) - (5|2|6)) - <5 56 ((5|1|2)(2|5|6) - (5|6|1)(1|2|6))) 



+ 1(2|(1 + 4)|3)(5|(1 + 2)|6)<5 56 



(1|(2 + 3)|4)(1|(3 + 4)|2)A 3 

^| m (5l2,S34,S56) 



+ 



2 (1|(2 + 3)|4)A 3 

[13](12)(3|(1 + 2)|6) 2 Lo(^f) l[13] 2 (23)(l|(2 + 3)|6) 2 Li(^f) 



[5 6] ti23 <1 3) <3|(1 + 2)|4> ai 2 2 [5 6] t 123 (l|(2 + 3)|4) s 2 3 

tl23 



[13] (l|(2 + 3)|6)(2|(l + 3)|6)Lo^J [13] (12) (1|(2 + 3)|6)(3|(1 + 2)|6) Mjg) 
[5 6]ti 23 <l|(2 + 3)|4) s 23 [5 6]ti 2 3<13)<l|(2 + 3)|4) s 23 

1 [64] 2 (42) 2 t 123 Ll (^t) [64] 2 (42)i 123 L °l 



- 2 



2(2 3)[5 6](1|(2 + 3)|4) t 2 23 [56] (1|(2 + 3)|4)(3|(1 + 2)|4) s 56 

[64] (42) (2|(l + 3)|6)Lo(^f) 
(23)[56](1|(2 + 3)|4) s 56 

1 \ ( _ / (12)[61][6 2]t 123 [12](51)(5 2)t 124 \ 

(3|(l + 2)|4)(l|(3 + 4)|2)A 3 L lSl2 S34j l [56] + (56) J 

+ (2si2-5 56 )(5 2) (1|(3 + 4)|2) [16] 

+ ((«13 + S23)(>23 + S 2 i) ~ S 12 (s 12 + S 23 - Si 4 ))(5|l|6) 

+ ((Sl4 + S24)(S13 + Sl4) - Sl 2 («12 ~ «23 + Sw)) <5 1 2 1 6> ln( — ^ ) 

J v— S56 7 

-flip, 



(3.15) 



where 'flip' is to be applied to all the preceding terms in F + . 
3.3 Subleading Color Primitive Amplitude 

Here we give the primitive amplitude Af(l, 2, 3, 4), which contributes only at subleading order 
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in N c . The "tree amplitude" appearing in eq. (3.3) is 



4 ree ' sl (l,2,3,4) 



[13] (54) <2| (1 + 3)|6) <24>[6 3]<5|(2 + 4)|1> 



■Sl2S56*123 



Sl2'S56*124 



(3.16) 



Note that 4 ree ' sl (2,3,l,4) = -^ ree ' ++ (l,2,3,4). 



tree, H — h / 



For the subleading-color primitive amplitude, it is convenient to introduce an 'exchange' op- 
eration where the 5,6 fermion pair is exchanged with the 1,2 fermion pair, 



exchange: 



1 



2^5. 



(3.17) 



The box-parent part of Af is even under this exchange. It is also convenient to define a symmetry 
operation 'flip_sl' (distinct from the leading-color 'flip'), 



flip_sl: 1 <-> 5 , 2^6, 3^4, (a b) <-> [a b] , (a | b |c ) <-> (a + | b |c + ) = (c | 6 |a ) . 

(3.18) 

The singular contribution is 



V sl (l,2,3,4) = 



^y_ 4 i + Li ^y_i /j^_y_7i 

2e v-W J L e 2 v-5i 2 y 2 e V-si 2 y 2j ' 

(3.19) 

where the first bracket is from the 'box parent' graphs in fig. 3a and the second bracket is from the 
'triangle parent' graphs in fig. 3b. In fact, the entire contribution from the triangle parent diagrams 
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is contained in the second bracket. The finite contribution is 



l2 / A r \ 2 

F sl (l,2,3,4) 



[13]" (45)" (3|(1 + 2)|6) 2 (2|(1 + 3)|4)^ ~2m/, 

J-iS_ 1 ^34, Ei23) «12) «56 J 



[1 2] (5 6) t 123 (4|(l + 2)|3) (1 2) [56] t 123 (3|(l + 2)|4)3 

+ T/| m ( Sl2 , S3 4,S56) + 



+ 



1 [6 4] 2 (4 2) 2 t 123 L i(^S) +2 [64](42)(2|(l + 3)|6)Lo(^f) 



2<12> [56] <3|(1 + 2)|4> i 2 23 (1 2) [56] <3|(1 + 2)|4) S56 



(2 3)(24)[6 4] 2 t 123 1 (2 3)[64](2|(1 + 3)|6) ( -t 123 \ 

(12)[5 6](3|(1 + 2)|4) 2 S56 2(12)[56](3|(1 + 2)|4) 2 U \ -s 5 J 



~, wT^r lnT— — — — ^ - flip.sl 
4 <1 2> [56] ti23<3|(l + 2)|4> V -s 56 J 

lu , ^2^3 [12] 



2 (3|(1 + 2)|4)AI ([5 6] (2 5) <534 ' 556 " 2 (1 2) (5 6) [6 1] 5l2 + 4 (1 2) S56 ( 5 2) [6 1] * 56 ) 
((5 2) (ai2-S34)-2<56) [61] (12)) 



-s 3 4 

1 [12] (2 5) 



2(3|(1 + 2)|4)(5 6) A 3 

I (3|(1 + 2)|4)A 2 ((5 2) [1 6] (534<556 + 2Sl2<5l2) + 2 ([1 2] [5 6] (2 5)2 + (1 2) (5 6) [6 1]2) ^ 
1 t 123 / (52)[16] _[12](2 5) 2 (12) [6 1] 2 



2(3|(1 + 2)|4)A 3 V 11 1 (5 6) [5 6] 

/ [12] [64] (3 2) (5|(2 + 3)|1) x 

+ l-(3|(l + 2)|4) 2 A 3 + (56) (3|(1 + 2)|4)A 3 J ((52) ~ 2 (5|(3 + 4)|1) (1 2)) 
[64] (3 2)t 123 



[56] (3|(1+2)|4) 2 A 3 



([16] <y 56 -2<2|(3 + 4)|6) [12]) 



+ (3|(1+ 4 2)|4)A 3 ((54) [41] (2|(1 + 3)|6) + [63] (23) (5|(2 + 4)1 1}) 
+ 2 (12) [56^(31(1 + 2)14)^3 ( " 2 (2|5|6)556 + (2 ' (3 + 4)|6) ' 34) 

" 2 (12) (56) (31(1 + 2)14)^ ^ « 53 > ^ ^ + < 54 > ^ ^ " W + « < 32 » 
[64] (2|(1 + 3)|6)(2"|(5 + 6)(1 + 2)|3+) + 1 (23) [46] (2|(1 + 3)|6) 



(12)[56]t 123 (3|(l + 2)|4) 2 2(12)[5 6](3|(1 + 2)|4) 2 

1 (2|(1 + 3)|6) 2 (52)((45)(2l(l + 3)|4)-(25)t 123 ) \ " 

4 (1 2) [5 6] t 123 (3| (1 + 2)|4) + (1 2) (5 6) i 123 (3| (1 + 2)|4) J P " 



1 (^123^34 + 2si 2 g 56 ) / [61] 2 (5 2)" x _ [61] (5 2) 

2 (3|(1 + 2)|4)A 3 Ml 2] [5 6] (1 2) (5 6)/ 1 123 124 } (3|(1 + 2)|4)A 3 



+ exchange, 

(3.20) 

where 'exchange' is to be applied to all the preceding terms in F sl , but 'fhp_sP is to be applied only 
to the terms inside the brackets ([ ]) in which it appears. The three-mass triangle coefficient T is 
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given by 



T-3s (t 6 + 2s s ^ ([12] (2 5) 2 [5 6] + (5 6) [61] 2 (12)) (f 128 - t 124 ) 

" 3 S3 4(ti23<J34 + 2s 12 s 56 ) ( 3 |(l + 2)|4)A§ 6 Sl2S56 s 3 4 (5 2) [6 1] (3|(1 + 9)|4)A 2 

■ ((5 2) [6 1] ( S56 + si2 + *34) + [1 2] (2 5) 2 [5 6] + (5 6) [6 l] 2 (1 2)) 



(3|(1 + 2)|4)A 3 



^ 2 + W-([661<5 2 >* B a-[61]<12>^) ^^^ 3 + ^ 



(3|(1 + 2)|4)A 3 U 1 ) (3|(1+2)|4)A 3 

( 3| ( f+2^ ] 2 A3 « 5 6 ) « 6 !] - t 6 5 1 < 5 2 ) P !]) " t 1 2 ] M<2 5) im - (2 1) [1 6] (6 5))) 
(23)[64](5|(2 + 3)|1) 1 (23) 2 [64] 2 

:(ll23C>34 + 2Si 2 S56j 



(3|(1 + 2)|4) 2 2 v " a M iZ ow (12)[5 6](3|(l + 2)|4) 3 

-2([6q<62)i M -[61]<12)* 12 )J^|- 

+ 2 (-(2|3|4)(4|5|6)<5 56 + <2|1|3)<3|4|6)<5 12 + (2|(1 + 3)|6) S34 <5 34 ) ' ''* )|(,) 



(12) [5 6] (3|(1 + 2)|4)A 3 



(2|(l + 3)|6)(5|(2 + 3)|l)s 3 4 ,. r e , (2|(1 + 3)|6) [14] (3 5) 

+ 4 /QlTl 1 oTTl\ A 1*123^34 + ^Sl2S56j /Q|/i I oM/l\2 

(3|(1 + 2)|4)A 3 si 2 s 56 (3|(l + 2)|4) i 

+ 2 (-[12] (2 3) [3 4] (4 5) + (5 6) [64] (4 3) [31] + (5|(2 + 3)|l)s 34 ) ^ ' + 

1. [16] (2 5) 

- 7r(*123C>34 + 2si 2 s 56 



S12S 5 6<3|(1 + 2)|4> 



2 v— . ■ ™^ 12S56 (3|(1 + 2)|4) 
~ I ^J^Kl + 2 2)|4) (2<2|(1 + 3)|6)«5|3|1) - <5|4|1» 

+ ((2|3|6) - (2|4|6))(5|(2 + 4)|1) + [61] (2 5)<5 34 ) 

+ ( /™ii\ + Kuii\+ 1 r, / ^ [14](35) i [13](45) 

+ f — <5|3|1> + (5|4|1) + -(t 123 - t^){ + m + m ) ) — — - 



(4(2|(l + 3)|6)-(2|3|6) + (2|4|6)). [1:!](l;,) 



S12S56 



(3.21) 



3.4 Axial Vector Contribution 



The primitive amplitude Ag x (l,2,3,4) is unique in that neither of the external quark pairs 
couples directly to the vector boson; instead they couple through the fermion-loop triangle diagram 
shown in fig. 4. The contribution to the amplitude (2.6) vanishes when the vector boson is a photon 
(by Furry's theorem, i.e. charge conjugation invariance). The Z contribution is proportional to 
the axial vector coupling of the Z to the quarks in the loop. As the u, d, s, c quark masses may 
be neglected, only the t, b quark pair survives an isodoublet cancellation in the loop, due to its 
large mass splitting. We use the results of ref. [31] to obtain this contribution; that paper presents 
the fully off-shell Zgg vertex, and we need only contract it with the three fermion currents. The 
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infrared- and ultraviolet-finite result is 



Ag*(l,2,3,4) = 



2i f(m t ; s 12 , s 3 4, g 56 ) ~ f(m b ; s 12 , s 3A , s 56 ) f [6 3] (4 2) (2 5) [6 1] [1 3] (4 5) 
"(4vr)2 S56 V (12) [12] 

+ (1 ~ 3, 2^4), 



) 



(3.22) 



where the integral f(m) is defined in the appendix. We need only the large mass expansion (for 
m = m t ) and the m = limit (for m = nib) of this integral; these are presented in the appendix. 



In this paper we presented all one-loop four-quark amplitudes which enter into the computation 
of the next-to-leading order QCD corrections to e + e~ — ► (j*,Z) — > 4 jets. The two-quark two- 
gluon amplitudes will be presented elsewhere [8,9]. We obtained the amplitudes by demanding that 
their functional forms satisfy unitarity and factorization. We also used a color decomposition and 
a helicity basis to express the amplitudes in a compact form. This way of obtaining amplitudes is 
significantly more efficient than using Feynman diagrams, since previously computed amplitudes 
are used to construct new ones. These amplitudes can be incorporated into numerical jet programs, 
which should lead to an improved knowledge of the QCD background to searches for new physics 
in various processes. Indeed, the leading-in-iV c parts of the contributions presented here, together 
with those for two quarks and two gluons [8,9] have already been inserted into one such program 
for four-jet production in e + e~ annihilation [10]. 
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We collect here the integral functions appearing in the text, which contain all logarithms and 
dilogarithms present in the amplitudes. Except for the contribution of the top quark to the fermion 
triangle loop in Aq x , all internal lines are taken to be massless. The following functions already 
appear in the evaluation of pentagon loop integrals where all external legs are massless, and of box 



4. Summary and Conclusions 



Appendix I. Integral Functions 
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integrals with one external mass, e.g. as occur in the one-loop five-gluon amplitudes [13]: 

t ( \ t < \ L o( r ) + 1 

L o(r) = -—, Li(r)= , 

1 r 1 r 2 (1.1) 

7T 

Ls_i(ri,r 2 ) = Li 2 (l - n) + Li 2 (l - r 2 ) + Inn hir 2 - — , 
where the dilogarithm is 

Ibto—r + Wzlt. (L2) 

Jo y 

The function Ls_i is simply related to the scalar box integral with one external mass, evaluated in 
six space-time dimensions where it is infrared- and ultraviolet-finite. 

The box function analogous to Ls_i, but for two adjacent external masses, is 

Ls 2 _f (s, t;mlml) = -uJl-^)-uJl-^)-W(^) + 1 -\n(^) In 



t J V t J 2 \-t J 2 \-mjJ \-m 



+ 



-(s — m\ — mj) + 



l ' ' - m?m 'j/| m ( S ,m?,m|), 

(1.3) 

where /| m is the three-mass scalar triangle integral. Here we use only a version of this box function 
with l| m removed, 

ST ( , „ ™i, = - 14 (1 - ^) - L , (1 - =1) - i * (^) + i ta (^) m (^) . 

(1-4) 

The analytic properties of these integrals are straightforward to obtain from the prescription 
of adding a small positive imaginary part to each invariant, — ► s^- + ie. One expands the 
logarithmic ratios, ln(r) = ln(^p-) = ln(— s) — ln(— s'), and then uses 

ln(— s — ie) = In \s\ — iirO(s) . (1.5) 

where O(s) is the step function: 0(s > 0) = 1 and G(s < 0) = 0. The imaginary part of the 
dilogarithm Li 2 (l — r) is given in terms of the logarithmic ratio, 

ImLi 2 (l — r) = — ln(l — r) Imln(r). (1.6) 

For r > the real part of Li 2 (l — r) is given directly by eq. (1.2). For r < one may use [33] 

IT 2 

ReLi 2 (l-r) = y - In \r\ In |1 - r\ - ReLi 2 (r) , (1.7) 



with ReLi 2 (r) given by eq. (1.2). 
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The analytic structure of I 3 m is more complicated [34,35,36], and the numerical representation 
we use depends on the kinematics. The integral is defined by 



^| m (si2,S34,S56) = / d 3 a,i5(l - ai - a 2 - a 3 ) 
Jo 



-s 12 aia 2 - s 34 a 2 a3 - s 56 a 3 ai 



(1.8) 



This integral is symmetric under any permutation of its three arguments, and acquires a minus 
sign when the signs of all three arguments are simultaneously reversed. Therefore we only have to 
consider two cases, 

1. The Euclidean region s\ 2 , «34,S56 < 0, which is related by the sign flip to the pure Minkowski 
region (si2, S34, S56 > 0) relevant for e + e~ annihilation. Here the imaginary part vanishes. This 
region has two subcases, depending on the sign of the Gram determinant A 3 (si2, S34, s^q) defined 
in eq. (3.6): 

la. A 3 < 0, 
lb. A 3 > 0. 

2. The mixed region s 12 , S56 < 0, s 34 > 0, for which A 3 is always positive. 

In region la one may use a symmetric representation found by Lu and Perez [34], which is 
closely related to that given in ref. [36]: 



r3m 



Cl 2 2tan" 



i/v^a; 



'12 



+ Cl 2 2 tan 



J34 



+ Cl 2 2 tan 



ifV=Al 



^56 



where the Sij are defined in eq. (3.6) and the Clausen function Cl 2 (x) is defined by 

sin(nx) 



C\ 2 (x) ee Yl 



n=l 



dt ln(|2sin(t/2)|). 



(1.9) 



(1.10) 



r3m 



In regions lb and 2 a convenient representation is given by Ussyukina and Davydychev [35], 
1 



Re 



2 (Li 2 (-px) + Li 2 (-py)) + ln(px) ln(py) + In (|) In ( - 



i ±py \ + t[_ 

+ px J 3 
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mt8(s34) f (612 + v / A 3 ")(<5 56 + 



(612 - V / A 3 ")(5. 



56 



As)/ ' 



(1.11) 



where 



2s 



56 



Sl2 S34 

x = — , y = — , p = 

S56 s 5e <5 56 + sj A 3 



(1.12) 



Finally, in the top quark contribution to Aq x the combination f(m t ) — f(rrib) appears, where 
f(m) is the integral 



/(m;si2, 534,555) = / d 3 ai5(l - ai - a 2 - 
Jo 



03) 



0203 



- si 2 aia 2 - s M a 2 a 3 - £560301 



(1.13) 
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This integral is complicated for arbitrary mass m; however, the large and small mass limits of it 
suffice for m t and m& respectively. For m = m t we simply Taylor expand the integrand in 1/m; for 
m = nib we set m& to zero, and reduce /(0) to a linear combination of the massless scalar triangle 
integral l| m given above, logarithms and rational functions. We get 

/K; , 12 , S3 4, .56) = 24^2 + ^4 + • • • . 

t(C\ ^ ( 3fi 34^34 1 \ r3 m/ \ . ( 3g 56^56 1 \ . f S 12 \ 

/(0; Sl2, S34, «56) = (^^2 ^"J Sl2S56^ 3 ( s 12, S34, S56J + 2A3" J V^S^/ 



/3Si2^12 1 \ , f-S5&\ #34 



(1.14) 
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